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Abstract 

<N 

t-H A fundamental limit to the application of Time Projection Chambers (TPCs) in high-rate experiments is 

the accumulation of slowly drifting ions in the active gas volume, which compromises the homogeneity of 
the drift field and hence the detector resolution. Conventionally, this problem is overcome by the use of 
Qh ion-gating structures. This method, however, introduces large dead times and restricts trigger rates to 
a few hundred per second. The ion gate can be eliminated from the setup by the use of Gas Electron 
Multiplier (GEM) foils for gas amplification, which intrinsically suppress the backflow of ions. This makes 

CO the continuous operation of a TPC at high rates feasible. 

In this work, Monte Carlo simulations of the buildup of ion space charge in a GEM-based TPC and the 
correction of the resulting drift distortions are discussed, based on realistic numbers for the ion backflow in 
a triple-GEM amplification stack. A TPC in the future panda experiment at FAIR, in which pp interaction 
rates up to 2-10 7 s _1 will be reached, serves as an example for the experimental environment. The simulations 
show that space charge densities up to 65 fC cm -3 are reached, leading to electron drift distortions of up 
to 10 mm. The application of a laser calibration system to correct these distortions is investigated. Based 
on full simulations of the detector physics and response, we show that it is possible to correct for the drift 

.^h distortions and to maintain the good momentum resolution of the GEM-TPC. 
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charge, drift distortions, laser calibration 
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1. Introduction 

A Time Projection Chamber (TPC) pQ can be 
regarded as an almost ideal device for charged- 
particle tracking. A large number of 3-dimensional 
hits measured along a particle track (typ. 50-100) 
eases the task of pattern recognition in a dense en- 
vironment and allows particle identification (PID) 
via the measurement of specific ionization. Large 
solid-angle coverage combined with very little ma- 
terial in the active part of the detector makes this 
device very attractive for applications in which high 
resolution is to be combined with small photon con- 
version probability and little multiple scattering. 
TPCs have been successfully used as large volume 
tracking devices in many particle physics experi- 
ments, e.g. PEP-4 [2], TOPAZ [3], DELPHI [4], 
ALEPH 0, NA49 ©, STAR [7], CERES [8] and 
ALICE [9]. 

In its standard form, a TPC consists of a large 
gas-filled cylindrical vessel sorrounding the inter- 
action point, placed inside a solenoidal magnetic 
field to measure the momentum of charged particles 
[TP] . The passage of an ionizing particle through 
the TPC produces a trace of electron-ion pairs. A 
uniform electric field along the cylinder axis ( "drift 
field"), created by the two endcaps and a field cage 
made of azimuthal metallic strips on the cylinder 
walls maintained at a linearly decreasing potential, 
separates positive gas ions and electrons. The ion- 
ization electrons then drift towards the readout an- 
ode located at one endcap of the cylinder. The 
transverse diffusion, which can become quite large 
for drift distances of the order of 1 m, is reduced by 
the magnetic field parallel to the drift direction. At 
the endcap, a plane of proportional wires is conven- 
tionally used for avalanche multiplication of ioniza- 
tion electrons. The signals induced on an arrange- 
ment of pad electrodes provide a measurement of 
the track projection onto the endplate. The third 
coordinate of the track is extracted from the mea- 
surement of the drift times of the ionization elec- 
trons. The reconstruction of a track from the mea- 
sured data requires a precise knowledge of the drift 
of electrons and hence of the electric and magnetic 
fields in the chamber. 

Electrons have drift velocities of the order of sev- 
eral cm |Lts _1 and are thus quickly removed from the 
drift volume. The drift velocity of ions, however, is 
three to four orders of magnitude smaller (see Table 
[l]for the values used for our simulations), leading 
to a slow buildup of space charge in the chamber. 



There are two principal sources of ions in a TPC: 

• Gas ionization by fast charged particles 
traversing the drift volume: The created ions 
slowly drift towards the cathode end-plate of 
the TPC. 

• Avalanche multiplication: During avalanche 
amplification, a large amount of electron-ion 
pairs are created, given by the total gain G, 
which is typically of the order of 10 3 to 10 4 . 
Without further measures, the ions created in 
the amplification process would move back into 
the drift volume and lead to significant distor- 
tions of the electric field. 

The total amount of charge accumulated in the drift 
volume depends on the rate and momentum distri- 
bution of the incident particles, the properties of 
the gas, and the amount of ions from the avalanche 
region drifting back into the drift volume. To pre- 
vent avalanche ions from reaching the drift volume, 
TPCs are normally operated in a pulsed mode, 
where an electrostatic gate to the readout region 
is opened only when an interaction in the target 
has occurred, and is closed immediately thereafter 
[TO] . The time needed to remove the ions as well 
as the switching time of the gate constitute dead 
times for the experiments, which limit the trigger 
rates to several hundred per second. 

Modern particle physics experiments, in contrast, 
require high interaction and trigger rates and little 
or no dead time of the detector systems. In order to 
benefit from the advantages of a TPC in a high-rate 
experiment, one has to find other means of space 
charge suppression. As an alternative to ion gat- 
ing, the usage of Gas Electron Multipliers (jGEMk ) 
[TT] for gas amplification has been proposed, since 
these devices feature an intrinsic suppression of the 
ion back-drift p~2j [13] . Operating a TPC at interac- 
tion rates which are large compared to the inverse 
drift time of electrons - which is of the order of 
100 |Lts for typical drift distances of 1 m - means 
that several events will be overlapping in one drift 
frame. The TPC hence acts as an "analog track 
pipeline" with signals arriving continuously at the 
readout pads. Instead of an event-based, triggered 
readout, the appropriate readout mode of such a 
device is then a continuous electronic sampling of 
signals combined with an autonomous detection of 
hits, which are further processed based on their in- 
dividual time stamps. The association of these hits 
to tracks and of tracks to distinct physics events 
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Figure 1: Geometry and dimensions of the active 
part of the panda TPC (cross section). 



( "event deconvolution" ) requires real-time tracking 
capabilities of the data acquisition system. The 
successful development of a continuously running 
TPC, though challenging, opens the possibility to 
benefit from the advantages of such a detector in 
future high-rate experiments. 

Among other options, a GEM-based TPC was 
proposed as the central tracker of the experiment 
panda (Antiproton Annihilations at Darmstadt) 
experiment [14] at the future international Facility 
for Antiproton and Ion Research (FAIR), panda 
will use an intense, cooled antiproton beam with 
momenta from 1.5 to 15GeV/c impinging on dif- 
ferent targets and will reach luminosities up to 
2 • 10 32 cm -2 s _1 , resulting in a pp interaction rate 
up to 2 • 10 T s _1 . The dimensions of the panda 
GEM-TPC are shown in Fig. [TJ It consists of 
cylindrical vessel with 150 cm length and an in- 
ner (outer) radius of 15.5 cm (41.5 cm) with GEM 
amplification at the upstream endcap [15], placed 
in a 2T solenoidal magnetic field. Owing to the 
fixed-target geometry of the experiment, the inter- 
action point is not located in the middle of cylinder 
axis, but shifted upstream. The TPC was shown 
to provide excellent standalone pattern recognition 
and momentum resolution of the order of a few 
percent [16]. The ungated, continuous operation 
mode of the TPC at the envisaged event rates at 
panda gives rise to about 4000 tracks which are 
superimposed in the drift volume at any given time. 

Track densities similar to this environment are 
expected for the ALICE TPC after an upgrade of 



the LHC, foreseen for the year 2018. At an expected 
luminosity for Pb-Pb collisions of 6 • 10 2T cm -2 s _1 , 
a continuous readout of the TPC is required in 
order to make full use of the interaction rate of 
about 50 • 10 3 s _1 . This cannot be achieved with 
the present gated MWPC-based amplification sys- 
tem. A replacement by a triple- GEM amplification 
using large-size foils is currently under evaluation. 

The ion leakage from a multi-GEM detector, 
however, is considerably larger than that from a 
closed gating grid, which is typically < 10 -4 [9]. It 
is therefore important to develop an understand- 
ing of the effects of residual space charge in a 
GEM-based, continuously running TPC. To this 
end, we have developed a computer simulation of 
space charge accumulation and drift distortions for 
a GEM-based TPC. The simulation is based on the 
following steps: 

• transport particles from minimum bias physics 
events through the detector setup and calcu- 
late their energy loss; 

• model the drift of ions to obtain the spatial 
distribution of space charge; 

• calculate the resulting electric field in the drift 
volume using finite element methods; 

• solve the drift equation for electrons in electric 
and magnetic fields to get a map of drift distor- 
tions as a function of the point of generation; 

• simulate a laser calibration system to measure 
the drift distortions. 

Section [2] describes the detector setup and de- 
fines the parameters used in the simulations. In 
Sec. [3] we describe how the charge accumulation is 
modeled to derive a space charge distribution in the 
drift volume of the TPC. Section 0] details how drift 
distortions in the inhomogeneous drift fields are cal- 
culated. Values for expected drift distortions in the 
simulated panda environment are given. Finally, 
in Sec. [5] we demonstrate how the drift distortions 
can be measured and corrected by means of an ar- 
ray of ionizing laser beams. 



2. Ion Backflow in a GEM-based TPC 

The Gas Electron Multiplier (GEM) [TQ consists 
of a 50 |Ltm thin insulating Polyimide foil with Cu- 
coated surfaces, typically 5 |Ltm thick. The foil is 
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perforated by photo-lithographic processing, form- 
ing a dense, regular pattern of (double-conical) 
holes. For standard GEM foils, the holes have an in- 
ner diameter of ~ 50 \im and a pitch of 140 \im. The 
small dimensions of the amplification structures 
lead to very large field strengths (9(50 kV cm -1 ) in- 
side the holes of the GEM foil when a moderate 
voltage difference of typically 300 - 400 V is applied 
between the metal layers, sufficient for avalanche 
creation inside the GEM holes. 

The dynamics of charge movement and avalanche 
creation inside the GEM holes are complicated. 
Figure [2] shows a Garfield / Magboltz [17] simula- 
tion qualitatively explaining the suppression of ion 
backflow from the amplification region. Two elec- 
trons (light lines) are guided into the GEM hole 
by the drift field (here: 250 V cm -1 ), and produce 
avalanches by ionizing gas molecules (dots). The 
ions created in the avalanches (dark lines) closely 
follow the electric field lines because of their much 
smaller diffusion. Most of the ions are collected on 
the top side of the GEM foil, because the field in- 
side the GEM hole is much higher than the field 
above the hole. Only a few ions drift back into the 
drift volume. The extraction of electrons from the 
hole is facilitated by a higher transfer field below 
the GEM (here: 3.75 kV cm -1 ). The electrons can 
then be transferred to another amplification stage 
or collected at the anode. 

Typically three GEM foils are combined in a 
stack, leading to effective gains of the order of 10 3 
- 10 4 and at the same time guaranteeing a stable 
operation without the occurrence of discharges [18] . 
The effective gain G e ff of the gas amplification sys- 
tem is given by the average number of electrons 
arriving at the readout anode divided by the 
number of ionization electrons Nr : 



G, 



eff 



(1) 



Gas ionization produces an equal amount of elec- 
trons and ions, hence Nf = = Nj. The effec- 
tive gain G e ff differs from the intrinsic gas gain of 
the amplification stage, as it already includes losses 
due to limited electron collection or extraction ef- 
ficiencies of the foils. These, as well as gain fluc- 
tuations, have not been explicitly simulated in the 
present work. 




Figure 2: Garfield / Magboltz simulation of charge 
dynamics of two arriving electrons in a GEM hole. 
Electron paths are shown as light lines, ion paths 
as dark lines. Spots mark places where ionization 
processes have occurred. The paths have been pro- 
jected on the cross section plane. 



We define the ion backflow a^3 



IB 



N7 



(2) 



where Nq denotes the number of positive ions ar- 
riving at the drift cathode. Note that this quantity 
also includes a contribution from ions created dur- 
ing the ionization process. Ion backflow values of 
IB = 0.25% have been reached experimentally in a 
magnetic field of 4 T in an Ar/CH 4 /C0 2 (93/5/2) 
mixture [T9] . 

A very convenient parameter to describe the ion 
backflow in the simulation is the ratio e of the num- 
ber of ions drifting back from the GEM amplifica- 
tion stage Nq and the number of ionization elec- 
trons Ni m . 



e = 



Ni 



(3) 



Since the number of ions arriving at the drift cath- 
ode is given by the sum of ions from ionization and 
the ones from the amplification stage, 



(4) 



1 We choose this quantity over other definitions in the 
literature since it can be easily measured via the ratio of 
cathode to anode current. 
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Table 1: Parameters of the space charge simulation. 



Event rate 


2 • 


10 7 s- 1 


Beam momentum 


Pp 


= 2.0GeV/c 


Gas mixture 


Ne/C0 2 (90/10) 


Average energy per 


Wi 


= 36.7eV 


ion pair 






Nominal drift field 


400 V cm" 1 


Magnetic field 


B = 


= (0,0,2.0T) 


Ion drift velocity 


u+ 


= 1.766 cm ms" 1 


Electron drift velocity 


u~ 


2.731cm M^s -1 


Longitudinal diffusion 


229 \xm cm" 1 / 2 


Transverse diffusion 


128 (xm cm -1 / 2 


TPC dimensions 


r = 


= 15.75. ..41.2 cm 


(active volume) 


z - 


= -39.5... 109.5 cm 


Ion backflow factor 


e — 


4 



the two quantities e and IB are linked by 
Ni + eTVi 1 + e 




IB 



G e ff 



(5) 
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Figure 3: Phase space distribution of primary par- 
ticles in the laboratory frame produced with the 
DPM generator for pp collisions at 2.0 GeV/c beam 
momentum. The plot shows the absolute momen- 
tum versus the cosine of the polar angle of 10, 000 
produced events. The sharp band of recoil protons 
from elastic scattering processes dominates the bulk 
of charged pions and kaons that are created in an- 
nihilation reactions and inelastic scattering. 



3. Simulation of Space Charge Buildup 



The suppression factor r\ defines the efficiency of 
ion-backflow suppression in the IGEMl stack: 



' N~ 7VlG eff G eff ' 

and thus 

e = ^G e ff . (7) 

As a rule of thumb, a suppression factor of n ~ 
1/Geff is usually considered necessary in order to 
operate a TPC without gating grid. 

The detector gas is a crucial parameter for the 
performance of a TPC [20] . In order to mini- 
mize multiple scattering and photon conversion, as 
well as space charge effects, a popular choice is a 
Ne/C02 (90/10) gas mixture [6 , 9 , which combines 
low diffusion even at moderate magnetic fields, low 
ionization density, large ion mobility and high ra- 
diation length. The results presented in this paper 
are based on realistic values for the ion backflow of 
IB = 0.25% and the effective gain of G eff = 2000, 
yielding e = 4 according to Eq. (|5|, i.e. four ions 
drifting back from the amplification stage per in- 
coming electron. The full set of simulation param- 
eters used for the presented results is summarized 
in Table[T] 



3.1. Ionization Charge Creation 

To model the spatial distribution of ionization 
charge in the TPC for the case of panda, the Dual 
Parton Model (DPM) [21 Monte Carlo event gener- 
ator is used to simulate antiproton-proton (pp) an- 
nihilations at an incident p momentum of 2 GeV/c. 
The resulting phase space distribution of charged 
particles is shown in Fig. [3] 

A large sample of these generated events is then 
passed through a GEANT [22] simulation of the 
detector geometry in order to calculate the energy 
loss of the particles in the TPC gas volume. For the 
simulations presented here we used GEANT3 with 
a simplified energy loss model as proposed in the 
ALICE TPC Technical Design Report [23 and ex- 
perimentally verified in [24]. From the GEANT en- 
ergy deposit AE(r, z,(j>) in a volume element dV 
around a point (r, z, <p) inside the TPC, summed 
up for a given number of simulated events N ev , 
the number of ionization electrons/ions Ni(r,z,(j)) 
produced in this volume element is calculated us- 
ing an average energy per electron-ion pair for 
the Ne/C0 2 (90/10) gas mixture of Wi = 36.7 eV 
[25] [26]. Here, r, z and (j) denote cylindrical coordi- 
nates, where the z axis coincides with the symme- 
try axis of the TPC. The creation rate of ionization 
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charge in dV is then 

Ni(r,z,<t>,t) = R(t) , (8) 

with R(t) being the event rate. 

3.2. Model of Space Charge Accumulation 

In general, the charge density p(r, z,(j),t) at time 
t in a volume element dV centered around a point 
(r, z, 0) inside the TPC volume is given by: 



p(r,z,</>,t) 



e 



(9) 



Here N(r, z, <fi, t') is the net rate of ions enter- 
ing/leaving the volume element dV in the time 
interval [£',£' + dt']. N(r, z,</>,tf) includes both 
ions directly created in gas ionization processes, 
Ni(r, z, 0, i), as well as ions moving through the vol- 
ume surface by drift and diffusion. 

Solving Eq. ([9| in principle requires perfect 
knowledge of the dynamics of charge creation and 
drift /diffusion inside the TPC volume. In order to 
arrive at a manageable, yet sufficiently realistic sim- 
ulation model, the following assumptions are made: 

1. Azimuthal symmetry: Owing to the cylindrical 
geometry of the TPC, we treat the problem 
in cylindrical symmetry so that the resulting 
charge density map can be represented in the 
(r, z)-plane. 

2. Constant luminosity: We assume that R(t) = 
const., i.e. the rate of charge created in gas ion- 
ization processes is constant on time scales of 
interest to us [O(10 us)]: Ni(r,z,t) = N\{r,z). 

3. Electrostatic forces between the ions are ne- 
glected: the ion drift proceeds along straight 
lines with constant velocity u + . 

4. The effect of diffusion on the motion of ions is 
neglected. 

In addition, we move from the infinitesimal vol- 
ume elements dV of Eq. ([9| to macroscopic vol- 
umes AV(ri,Zj) = AV lJ ' with constant bin widths 
Az = 1.0 cm and Ar = 0.98 cm. Each such bin (ij) 
then represents a ring-shaped volume in the TPC: 



AV l - 



(10) 



where r^ut = r* + |Ar and r ijin = 
outer and inner radius of bin (i,j), respectively. The 



\Ar are the 



integral in Eq. ([9| can then be solved numerically 
by a Riemann sum: 



P(r,z,t) 



N ijk At 



k=i 



AV ij 



(ii) 

with the time t given by the number of steps in time 
bins t = n • At. 

During each time step At, the ionization charge 
iVj J At is newly created in bin of the chamber. 
Now the backflow of ions from the amplification 
stage has to be included. As can be seen from Ta- 
ble [l] the drift velocity of electrons exceeds that of 
ions by several orders of magnitude. In order to cal- 
culate the ion space charge we can therefore assume 
instantaneous electron drift. As a consequence, a 
contribution proportional to the total amount of 
newly added gas ionization charge in the full TPC 
volume has to be added to the first bin in z (j = 1) 
for each time bin At. The total number of ions 
added in each time step is therefore 



N{ j = N{ j At + S ld eJ2 K 3 ' At 



(12) 



j' 



with Si j being the Kronecker delta symbol. 

As a last ingredient, the drift of ions has to 
be taken into account. For the geometry of the 
panda TPC, where the readout anode is located at 
the upstream endcap of the chamber, the ion drift 
proceeds along the positive z direction towards the 
cathode endcap. We choose At to correspond to 
the time needed for an ion of drift velocity u + to 
travel the distance Az. This allows us to include 
the ion drift in a time bin At simply by shifting the 
total charge by one bin in z. We therefore arrive at 
a recursive algorithm to calculate the space charge 
in the n th time step: 



N ijn = N ij + N i u 



(13) 



From the equations above it is clear that an equi- 
librium space charge in time, N ljn = 7V* J ' n_1 , will 
be reached everywhere in the chamber as soon as 
the perturbation due to ion backflow reaches the 
cathode, i.e. for t = 2; max /'U + , where z max is the 
maximum drift distance. Figure [4] shows the space 
charge distributions calculated with this method 
after the first step in time (i.e. immediately after 
switching on the beam) and after reaching equi- 
librium. The input parameters to the simulation 
are given in Table [l] After the first step, Fig. [4a, 
the distribution of ions from direct gas ionization 
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Figure 4: Space charge density p(z,r) in the 
panda TPC for 2 • 10 7 s _1 pp annihilations at 
beam momentum Pp = 2.0GeV/c with ion back- 
flow factor e = 4, using the recursive algorithm of 
Eq. (13). (a) After one step the distribution of 



ion charge from gas ionization processes is visible 
superimposed with the corresponding charge back- 
drifting from the amplification stage (red band at 
z = —40 cm). The band of high primary ioniza- 
tion density at z ~ is caused by slow protons 
from hadronic elastic pp-scattering. (b) Final equi- 
librium space charge density. 



by charged particles from pp annihilations with a 
band at a polar angle close to 90° from elastic 
pp-scattering is visible, together with the sheet of 
charge at z ~ —40 cm starting to drift back from the 
amplification region for this time bin. Due to the 
fixed target geometry of panda and the resulting 
forward boost of the scattered particles, the charge 
distribution is not symmetric in z around the inter- 



action point. Figure 4b shows the equilibrium space 
charge distribution which is reached after ions from 
the amplification stage have reached the cathode. 
The final space charge density varies between 20 
and 65 fC cm -3 . 



4. Drift Distortions 

4.1. Electric Field due to Space Charge 

From the charge-density distribution the electro- 
static potential ip is computed by solving Poisson's 
equation, which in cylindrical coordinates reads: 



A(p 



ld_ 
r dr 
1 



dip 
dr 

P(r, </>, z) 



1 d 2 p 



d 2 cp 
dz 2 



(14) 



In case of rotational symmetry, p(r, 0, z) = 
p(r,z), Eq. (14) reduces to 



dr 2 



I dip 
r dr 



d 2 y 

dz 2 



eo 



-p(r,z) 



(15) 



Equation ( 15 ) is solved numerically with the Finite 
Element Method. The corresponding variational 
formulation [27] is given by: 



/ 



Vp-Vv- 



ldip 
r or 



dzdr = 



h 



-dz dr 



(16) 



Here, ft is the drift region in the r-z-plane, V = 
(d/dr,d/dz), and v is an arbitrary, dimension- 
less test function which satisfies the same Dirichlet 
boundary condition as ip: 



(17) 



where dQ represents the boundary of the drift re- 
gion, i.e. the field cage wall. This boundary con- 
dition corresponds to a perfectly conducting field 
cage, such that the potential of the field-cage is 
fixed to zero. The software package Dolfin/FeniCs 
[28] is used to realize the finite element solver. 

By calculating the gradient of the resulting po- 
tential we finally obtain the electric field E caused 
by the space charge. Due to the assumed ip- 
symmetry, the non- vanishing components of E are 
the r and z components. The drift field = 
400 V cm -1 is oriented along the z direction and is 
superimposed onto the distortion field (see Fig. [5|. 



4.2. Electron Drift Offsets 

To compute a map of the drift distortions as a 
function of the starting point of the drift, the equa- 
tion of motion for electrons drifting in electric and 
magnetic fields (E, B), 



d 

m — u 

dt 



eE + e[u" x B] - Ku~ 



(18) 
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Figure 5: Electric field caused by accumulated 
space charge in the panda TPC for 2 • 10 T s _1 pp 
annihilations at beam momentum pp = 2.0GeV/c 
with ion back-flow factor e = 4. (a) Radial com- 
ponent E ri and (b) component along drift direc- 
tion E z superimposed with the nominal drift field 
E d = 400 V cm" 1 . 



is integrated with a fourth order Runge Kutta al- 
gorithm^ In Eq. ^S) m is the electron mass, u 
is the macroscopic electron drift velocity, e is the 
electron charge and K = ej ' pT is a constant friction 
term modeling the microscopic stop-and-go motion 
of the drifting electron due to collisions with the 
gas components (mobility pT = u~ /E z with E z 
being the z component of E). The electric field E 
is the result of the steps described in Sec. |4.1[ while 
the magnetic field B is given by the calculated field 
map of the panda solenoid, with a nominal field 
of 2 T and a relative uniformity of better than 2% 
everywhere in the active volume of the TPC [29] . 
The fact that both E and B have non-zero compo- 
nents perpendicular to z leads to a displacement of 
drifting electrons compared to the ideal drift with 
constant velocity along z. Figure [6] shows the ra- 
dial and azimuthal components of this deviation at 
the readout plane as a function of the coordinates 
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Drift distortions £ of electrons due to 
charge accumulation, (a) in radial direc- 
(b) in azimuthal direction. The graphs 
deviations from a straight line drift ex- 
by an electron which starts its drift at 
z) in the TPC. The z-component of this 
map is found to be negligible in compar- 



of the starting point of the drifting electron. It 
should be noted that although both fields are as- 
sumed to exhibit cylindrical symmetry, there are 
nevertheless drift distortions perpendicular to the 
r-z-plane. This is due to the E x B term in the 



solution to Eq. (18). 



2 When reading field values from maps with finite bin size, 
linear interpolation in between bin centers is applied. 



5. Measurement and Correction of Drift 
Distortions 

In order to correct track data for space charge ef- 
fects, the resulting drift distortions have to be mea- 
sured during the operation of the TPC. This can 
be achieved by creating a well defined pattern of 
electron sources in the TPC and analyze its image 
as measured by the detector. A comparison with 
the expected image directly yields the distortions. 
Straight line ionization tracks from UV lasers can 
provide such a pattern [30] and have already been 
used for calibration in other drift chambers (e.g. in 
STAR [31 ). To assess the potential of this method 
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in the case of an ungated TPC, we have simulated 
such a laser system. 

5.1. Laser Calibration System Simulation 

Having found the drift distortions in z-direction 
to be small, we want to measure the distortions in 
the r-(/>-plane as precisely as possible. Ideally, this 
can be done with laser tracks parallel to the drift 
direction / detector symmetry axis. We chose a 
grid of 6 laser beams aligned equidistantly in radial 
direction, where the innermost (outermost) laser 
is placed 0.5 cm away from the inner (outer) field 
cage wall. One such row of lasers is created ev- 
ery 20° in azimuthal direction. The laser beams 
are parametrized as straight line tracks with an av- 
erage ionization density of 45e~cm _1 [31], and a 
Gaussian beam profile (a — 300 |Ltm). 

In the simulation routines, the electrons created 
by the laser rays are then drifted to the readout 
anode, applying the calculated distortion map and 
taking into account diffusion. The induced signals 
are calculated using GEM amplification and sig- 
nal induction with a realistic response of hexagonal 
readout pads with an outer radius of 1.5 mm as de- 
signed for the panda TPC. In a last step, signals 
adjacent in space and time are combined to hits by 
a clustering algorithm. A sample of such hit data as 
written out by the simulation framework is shown 
in Fig. [7| 

5.2. Reconstruction of Drift Distortions 

The reconstruction of the laser tracks is greatly 
simplified by the knowledge of the geometry (and 
the time, at which the laser grid has been illumi- 
nated): The hits are assigned to tracks by apply- 
ing simple residual cuts, such that all hits inside 
a tube of reasonable radius (15 mm for the results 
presented) around the nominal track position are 
associated to the corresponding laser track. It is 
assured that each reconstructed hit is assigned to 
only one (the closest) track. For each hit, a residual 
vector is obtained in the plane perpendicular to the 
corresponding track. 

In order to map out drift distortions as a func- 
tion of the point of electron creation in the drift 
volume, the large amount of signal data from the 
laser calibration needs to be fitted, smoothed and 
parametrized. For this purpose we have imple- 
mented an algorithm for fitting of a bi-cubic spline 
surface s(z,r) [32 to the measured residual data 
(c.f. Fig. [7]). A two-dimensional mesh of h x k 
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Figure 7: Positions of TPC hits (after running of 
the clustering algorithm) from a laser event. The 
laser tracks are sketched as lines along the drift di- 
rection in the TPC. For better visibility, distortions 
have been scaled by a factor of 2 for this image. 
Note that only the radial component of the distor- 
tions can be seen in this projection. 

knots (Xi,fjLj) over the data area (z,r) is chosen F] 
At each of these knots two cubic B-splines M^(z), 
Nj(r) are attached, each of which has a fixed shape 
and spans over five mesh knots A^_4 5 ...^ (/ij-4, j). 
Coefficients 7^ at every knot scale the B-splines: 

/i+4 /e+4 

1=1 j=i 

Assuming the underlying noise (mainly given by 
diffusion during drift) to be Gaussian, the surface 
fitting problem is equivalent to finding the set of 
coefficients 7^ which minimize the sum of squared 
distances 

n r n r 

S = Y,[s(Zr,r r )-f r ] 2 = Y,R 2 r (20) 
r=l r=l 

to the residual data f r {z r ,r r ) (where r = l...n r ). 
The corresponding matrix-equation ( "normal equa- 
tion") reads: 



A T A 7 = A T f , 



(21) 



where A is an n r x (/i+4)(/c+4)-dimensional matrix 
containing the spline information, 7 is the vector of 



3 In addition 4 knots on each side outside the data area 
are required to define the full set of B-splines. 
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knot coefficients with (ft + 4)(fc + 4) elements and f 
is the vector containing the n r measured residuals. 
Weighting of hits can be incorporated by introduc- 
ing a weight matrix W. Under the assumption of 
independent hit errors, W becomes diagonal, and 



instead of Eq. ( 20 ) one now has to minimize 



Equation (21) then reads 



A' T A' 7 = A' T f , 



(22) 



(23) 



where A' = wA, f = wf , and w rr = \/W rr . The 
weights for each hit are obtained from the covari- 
ance ellipsoid available after clustering, taking into 
account the orientation of the laser track. By set- 
ting the error to be large along the track direction, 
it is possible to work also with laser grids that are 
not optimized to the distortions topology. 

The solution to Eq. ( [23] ) can be obtained 
by matrix inversion or, more stably, by QR- 
decomposition. Figure [8] shows the least squares 
spline solution to the reconstructed residual data 
for both radial and azimuthal drift distortions. One 
laser event as described in Sec. 15.11 is used, result- 



ing in ~ 35, 000 reconstructed residuals. Typical fit 
times are of the order of 1 s on a single CPU desktop 
PC. 

5.3. Quality of Reconstructed Distortion Map 

The choice of the number of spline-mesh knots 
h x k over the data area is very important for the 
procedure. On the one hand, a higher number of 
knots means higher overall fit precision, but over- 
fitting becomes an issue. A lower number of knots, 
on the other hand, results in a smoother result, but 
the fit might not be able to follow the important 
features of the data set. We found (h x k) = 5 x 3 
to be a balanced value. 

To be able to judge the quality of the result shown 
in Fig. [8| we compare the reconstructed distortion 
maps to our original drift distortion input (Fig. |6| 
in Fig. [9j which shows the distribution of differ- 
ences between the reconstructed and the calculated 
distortion maps. 

Fitting a Gaussian to the peak of the distribu- 
tions gives cTQauss ^ 50 |Ltm and a mean well cen- 
tered around 0, showing that the total systematic 
error of the laser correction method is small and 
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Figure 8: Reconstructed drift distortions £ in (a) 
radial direction and (b) azimuthal direction, based 
on one laser event (compare to Figs. 6a and 6b). 



the overall precision of the presented distortion re- 
construction and fitting method is better than the 
expected single-hit resolution of the TPC. Hence, 
using a grid of laser beams, drift distortions in the 
TPC can be measured with a precision better than 
(9(100 M-m), not taking into account mechanical dis- 
tortions of the calibration system. 

5.4- Effect on Track Reconstruction 

The most important question to be answered is of 
course the impact of space charge effects on track 
reconstruction. As explained in the previous sec- 
tions, the steps required to obtain a model of space 
charge and study its possible effects and their cor- 
rection are: 

1. calculate space charge based on background 
events; 

2. calculate the drift distortion map; 

3. create one "laser event" consisting of the full 
laser beam mesh; 

4. reconstruct the laser tracks and fit the data 
with a spline; 

5. use this spline to correct for drift distortions. 

In the running experiment, measurements of drift 
distortions have to be constantly updated, requiring 
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(a) 




<p, reco <p 



(b) 



Figure 9: Direct comparison of drift distortions as 
reconstructed with the Laser grid (spline fit) and 
original simulation input. Shown are the differences 
of the original (binned) deviation map and the fit 
evaluated at the bin center for the (a) radial and 
(b) azimuthal drift distortions. The corresponding 
r.m.s. of the histograms are 72 \im for (a) and 86 |Ltm 
for (b), respectively. 



steps 4 and 5 to be triggered on a timescale short 
enough to be able to resolve possible space charge 
fluctuations. 

To study the effect of drift distortions on charged 
particle tracks, we simulated a test sample of 5000 
pion tracks (7r + ) at a momentum of 0.5 GeV/c and 
uniformly distributed scattering angle (20° < 6 < 
115°) in the panda TPC. Tracks are reconstructed 
using a Kalman-filter-based track fitting framework 
[33 . Figure 10 visualizes the impact of uncorrected 
drift distortions on the momentum measurement 



compared to the ideal situation of a completely ho- 
mogeneous electrical field. Without correction, the 
distribution of reconstructed momenta is deformed 
and systematically shifted towards larger values by 
space charge effects. This can be understood within 
the scope of the preceding sections: Depending on 
its polar angle, a track is more or less affected by 
drift distortions (c.f. Fig. [6]), leading to the broad- 
ening of the distribution. The asymmetric nature of 
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Figure 10: Effect of drift distortions (and their 
correction) for 7r + -tracks at 0.5GeV/c momentum 
(primary and secondary GEANT tracks). Right 
peak: Uncorrected case. Left peak: Hit posi- 
tions corrected with spline fir results. Transparent 
distribution: Simulation with no drift distortions 
present. 



the drift distortions with respect to the radial coor- 
dinate causes the measured curvature of the track to 
appear smaller than it actually was, leading to the 
shift to higher momenta. Applying the correction 
derived from the laser system essentially recovers 
the momentum reconstruction, with a resolution of 
(Jp/p = 1.74%. For comparison, the distribution 
of reconstructed momenta for the ideal case with- 
out any space charge effects yielding a resolution of 
(Tp/p — 1.59% is shown as outline in Fig. [Toj The 
reconstructed mean changes from 5.01 • 10~^GeV/c 
in the ideal case to 4.99 • 10 _1 GeV/c, respectively. 



6. Conclusions and Outlook 

In this paper we have presented detailed studies 
on space-charge buildup in a continuously operating 
GEM-TPC without a gating grid exposed to high 
interaction rates. We have presented a recursive 
algorithm to calculate the space charge accumulat- 
ing in the drift region, starting from minimum bias 
physics events, and taking into account the primary 
ionization as well as the ions drifting back from 
the amplification region. Assuming constant lumi- 
nosity, an equilibrium space-charge distribution is 
reached after one full drift time of ions through the 
chamber. In the example of the environment of 
the future panda detector, ion space charges of up 
to 65 fC cm -3 are reached. The electric field re- 
sulting from this space charge is calculated using 
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a finite element method to solve the Poisson equa- 
tion. The drift of electrons in the full electric and 
magnetic field of the setup is calculated from the 
Langevin equation solved by a fourth order Runge- 
Kutta method. Drift-path distortions of (9(1 cm) 
are found for the PANDA enviromment, requiring 
corrections to be applied to the measured hit coor- 
dinates. 

To this end, we have simulated a laser calibra- 
tion system creating a regular grid of straight ion- 
ization tracks at known positions in the chamber. 
The geometry of the laser grid was optimized to 
yield an accurate measurement of the expected dis- 
tortions, which are dominant ly in radial and az- 
imuthal direction for the panda case. In the case 
of azimuthal symmetry of the distortions, a two- 
dimensional spline fit to the measured residuals is 
sufficiently fast and accurate, reproducing distor- 
tions with a precision better than 100 |Ltm, thus 
opening the possibility to be applied even in the on- 
line reconstruction software. This distortion map 
is then used to correct hits from charged particle 
tracks. Using a Kalman-Filter-based track fitting 
we are able to reconstruct the momentum of tracks 
with a resolution very close the ideal case of no dis- 
tortions in the chamber. 

The method to calculate the space charge effects 
and their correction presented in this paper is also 
applicable to the case of other TPCs with azimuthal 
symmetry, e.g. at STAR or ALICE. If one cannot 
exploit symmetry features of the distortions or it is 
not possible to realize the optimal laser grid geom- 
etry due to technical constraints, proper treatment 
of errors taking into account the grid geometry, as 
realized in the present work, is important to mini- 
mize systematic errors arising from residual projec- 
tions. This allows the algorithm to be extended to 
make use of track information from other detectors 
instead of a fixed laser grid to measure the drift 
distortions. 
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